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Abstract 

Non-leptonic kaon decays are often described through an effective chiral weak Hamiltonian, 
whose couplings ("low-energy constants") encode all non-perturbative QCD physics. It has 
recently been suggested that these low-energy constants could be determined at finite volumes 
by matching the non-perturbatively measured three-point correlation functions between the 
weak Hamiltonian and two left-handed flavour currents, to analytic predictions following from 
chiral perturbation theory. Here we complete the analytic side in two respects: by inspecting 
how small ( "e-regime" ) and intermediate or large ( "p-regime" ) quark masses connect to each 
other, and by including in the discussion the two leading AI = 1/2 operators. We show 
that the e-regime offers a straightforward strategy for disentangling the coefficients of the 
AI =1/2 operators, and that in the p-regime finite- volume effects are significant in these 
observables once the pseudoscalar mass M and the box length L are in the regime ML < 5.0. 
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1. Introduction 



Understanding why the AI = 1/2 amplitudes for non-leptonic kaon decays are so much larger 
than the AI = 3/2 amplitudes, is a long-standing problem for QCD phenomenology. It has 
been known since the early 70s that the bulk of the enhancement must be due to strong 
interactions at low energies [1]. Therefore a reliable explanation must eventually be based on 
systematic non-perturbative methods, in particular on lattice QCD [2, 3]. 

It was realized long ago that instead of computing directly the decay amplitudes with 
lattice QCD, a simpler alternative is to use lattice simulations to determine the relevant low- 
energy constants (LECs) of the effective chiral weak Hamiltonian that describes kaon decays 
[3], and then use chiral perturbation theory to compute the physical amplitudes [3]-[8]. The 
determination of the LECs can be achieved by matching certain observables computed in 
lattice QCD and in chiral perturbation theory (%PT), as close as possible to the chiral limit. 
In this respect it is advantageous to approach the chiral limit by first extrapolating to small 
or zero quark masses, and increase the volume only afterwards. This setup corresponds to the 
so-called e-regime of xPT [9] (see also Ref. [10]). The power-counting rules in this regime [9] 
guarantee that the contamination from higher order LECs is reduced very significantly. In 
other words, the number of LECs that appear at the next-to-leading order (NLO) in the e- 
regime of xPT is typically much smaller than that at the next-to-leading order in the standard 
p-regime, where the infrared cutoff is provided by the pion mass rather than the volume. 

The matching of lattice QCD and the chiral effective theory in the e-regime has recently 
been considered in order to extract the strong interaction LECs [11]— [18]. Subsequently, it has 
been pursued for the determination of the weak LECs that we are interested in [15, 19, 20], 
as well as for the study of baryon properties [21]. This progress has been possible thanks 
to the advent of Ginsparg- Wilson formulations of lattice fermions [22]-[29], which possess an 
exact chiral symmetry in the limit of vanishing quark masses. Simulations in this regime are 
however challenging on the numerical side, and Refs. [15, 17] introduced several important 
technical advances in order to make them possible. 

In Ref. [20], a strategy based on these methods has been proposed to reveal the role that 
the charm quark mass plays in the AI = 1/2 rule. In particular, following the suggestion 
of Ref. [15], the observables that are considered are three-point correlation functions of two 
left-handed flavour currents and the weak operators. The first step is the matching of these 
observables, to extract the LECs of the weak chiral effective Hamiltonian, in a theory with 
a light charm quark, that is in a four-flavour theory with an exact SU(4) symmetry in the 
valence sector. The results of this computation can be found in Ref. [30]. The next step of 
the strategy is to increase the charm quark mass and monitor the LECs as we move towards 
a theory with an SU(3) flavour symmetry [20, 31]. 

In a previous paper [19], we have already computed the NLO e-regime predictions for the 
correlators of left-handed flavour currents and the AI = 3/2 weak operator, whose coefficient 
determines the kaon mixing parameter Bk in the chiral limit. The purpose of the present 
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paper is to extend the results of Ref. [19] in two ways. First of all, we compute the same 
observables as before, but also at larger quark masses, corresponding to the p-regime of chiral 
perturbation theory. The goal is to obtain a better understanding of the regions of validity 
of the e and p-regimes. Second, we include the AI = 1/2 weak operators in the analysis. 

We find that the e-regime does offer a clean way of disentangling the coefficients of the two 
leading-order AI = 1/2 operators. 

It is well known that the description of quenched simulations, which still are widely in use 
today, through a quenched version of chiral perturbation theory, is rather problematic. In 
particular the p-regime is strongly affected by quenched ambiguities that increase significantly 
the number of LECs [32], making it difficult to identify those that should be closest to the 
ones in the full theory. We have studied the effect of these ambiguities also in the e-regime 
at NLO, and find that they are significantly less severe in this case. 

In most of our analysis we will concentrate, however, on the full physical theory. The most 
immediate applications might then follow through the use of mixed fermion frameworks [33] , 
though progress towards dynamical Ginsparg- Wilson fermions is also taking place [34]. 

It should be made clear from the onset that choosing to consider correlators involving left- 
handed flavour currents in this paper, is not meant to indicate that they would necessarily 
be the ultimate way for determining the weak LECs. For instance, employing the zero-mode 
wave functions of the massless Dirac operator might also lead to a useful probe, even though 
for the pion decay constant they seem to be slightly disfavoured in comparison with the 
left-handed flavour currents [16]. 

Other methods to obtain the weak LECs have also been considered in the literature. For 
lattice approaches without an exact chiral invariance see, e.g., the recent work in Refs. [35]. 
For models inspired by the large- iV c expansion see, e.g., Refs. [36, 37]. 

This paper is organised as follows. We formulate the problem in Sec. 2, discuss the various 
regimes of chiral perturbation theory in Sec. 3, address the AI = 3/2 operators in Sec. 4, 
and the AI = 1/2 operators in Sec. 5. We conclude in Sec. 6. 

2. Formulation of the problem 

We start by considering QCD with 4 flavours. The quark part of the Euclidean continuum 
Lagrangian reads 

4 

Le = ^tprilnDn + m r )tp r , (2.1) 

r=l 

where r is a flavour index; the Dirac matrices 7^ are assumed normalised such that 7^ = 7^ , 
{7^5 lu] = 2$iiu', is the covariant derivative; m r is the quark mass; colour and spinor indices 
are assumed contracted; and repeated indices are summed over, even when no summation 
symbol is shown explicitly. In the following we will consider the three lightest quarks as 
degenerate in mass, m u = = m s = m, while the charm quark is heavier, m c S> m. 



2 



After an operator product expansion in the inverse W boson mass, weak interactions can 
be described with the Fermi theory involving four-quark operators. In the CP conserving case 
of two generations, the effective weak Hamiltonian is then [1] (for reviews see, e.g., [38, 39]) 

H w = 2V2G F V vd V: s \j2 hl([O w ]: uud - [O w Y ccd ] + h m [O m ] sd \ +R.c. , (2.2) 

where h^,h m are scheme-dependent dimensionless Wilson coefficients, with leading order 
values = l,h m = 0. The coefficients /i^ are known to two loops in perturbation theory 
[40], while h m remains undetermined. In Eq. (2.2) we have introduced the notation 

(2.3) 

P w ] rsuv = (^P-M^P-^v) , (2-4) 
[Orn] sd = (ml-m 2 u ){m s ^ s P^ d ) + m d ^ s P + ^ d )}. (2.5) 

Here r,s,u,v are generic flavour indices, while u,d,s,c denote the physical flavours. The 
chiral projection operators P± read P± = (1 ± 7s)/2, where 75 = 70717273- The colour and 
spinor indices are assumed to be contracted within the parentheses. 

In order to match the Hamiltonian of Eq. (2.2) to the one in the SU(3) chiral theory, the 
first step is to decompose it into irreducible representations of the SU(3)lxSU(3)r flavour 
group, present at low energies. The weak operators are singlets under SU(3)#, and projecting 
them onto irreducible representations of SU(3)l, the weak Hamiltonian can be rewritten as 

H w = 2v^G F V^{/i+[^ 

- ^(^ + ^)[O ro ]w-^(^-^)[O w ] SC(ic + /im[Om] s d|+H.c. , (2.6) 

where 

[dw]tuud = 2 {[^ w ] sW + [O w ]sudu - jr ^2 (\Qw]skdk + [O w ]skkd^ } , (2-7) 

k=u,d,s 

[Pw]fd = 2 X/ (\Qw]skdk ± [O w } s kkdj ■ (2.8) 
k=u,d,s 

The first operator in Eq. (2.6) transforms under the 27-plet of the SU(3)l subgroup: it is 
symmetric under the interchange of quark or antiquark indices, and traceless. The remaining 
ones, transforming as 3* (g> 3 and being traceless, belong to irreducible representations of 
dimension 8. 

If, as the next step, the charm quark is also integrated out, then the operators in Eq. (2.6) 
go over into the standard ones, commonly denoted by Qi, i = 1,...,6 [41, 42] (of which five 
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are independent). It is probably safer to keep the charm quark in the simulations, though, 
since integrating it out perturbatively is not guaranteed to be a safe procedure. Moreover, 
the quenched three-flavour theory contains spurious operators [32]. For these reasons, we 
prefer to consider the four-flavour theory of Eq. (2.6) to be the QCD-side of our problem. 

Now, at large distances, the physics of QCD can be reproduced by chiral perturbation 
theory. For a degenerate quark mass matrix, the leading order chiral Lagrangian reads 



^ X PT = — Tr 



F ~ - 'd.Ud^] - ^Tr[e ie ^U + U^'^] , (2.9) 



where U £ SU(iVf), iVf = 3, and 9 is the vacuum angle. Apart from 9, this Lagrangian 
contains two parameters, the pseudoscalar decay constant Fand the chiral condensate S. At 
the next-to-leading order in the momentum expansion, additional operators appear in the 
chiral Lagrangian, with the associated low-energy constants Li,L 2 , ... [43]. 

Obviously the chiral model can be extended to include a weak Hamiltonian [44]. We denote 
the chiral analogue of H w in Eq. (2.6) by TL W . To again define dimensionless coefficients, we 
write TL W in the form [3, 5] 

H w = 2V2G F V ud V: s ^g 2 70 27 + 2g 8 8 + 2g' 8 0' 8 \ + H.c. , (2.10) 

where 527,58 and d's are the low-energy constants we are interested in. The operators read 

027 = [6w\ + suud = \ ) (pw\ sudu + ^[°w] S uud) > t 2 - 11 ) 

P w ] rsuv = tj-faurf) (dJJlA , (2.12) 

4 V / ur V / vs 

Os = [K w ]f d = l Yl [°»U M , (2-13) 

k=u,d,s 

> 8 = ^mzfeWKU + rfe-VM) , (2.14) 

2 V / ds 

where we have made use of Tr [d^UU^] = to simplify the chiral versions of Eqs. (2.7), (2.8). 

In the following, we will find it useful to generalize the notation somewhat from the standard 
SU(3) case introduced above. Let N v = 3 be the number of valence flavours, and N{ the 
number of degenerate sea flavours in the chiral Lagrangian. The standard case corresponds 
to Nf = N v , but one can also envisage other interesting situations, for instance Nf = 4 [20], 
or Nf — > 0. We note that the simplified forms in Eqs. (2.11), (2.13) only apply for iVf = A^ v ; 
in general, the combinations in Eqs. (2.7), (2.8) need to be employed (the generalizations of 
these combinations to arbitrary 7V V , Nf are summarised in Appendix A). In the remainder of 
this Section we have in mind the case Nf = N v but the formulae are written in a way which 
will be useful in Appendix C, where we analyse the situation Nf / iV v . 

The principal strategy now is to construct three-point functions by correlating H w with 
two left-handed flavour currents on the QCD side, and to match to predictions from xPT for 
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the same objects. In QCD, the left-handed flavour current can formally be defined as 

Jl ee #T a 7/i P_^ , (2.15) 

where T a is a traceless generator of the valence group SU(iV v ), and all colour, flavour, and 
spinor indices are assumed contracted. Note that J® defined this way is formally purely 
imaginary. 3 

The two and three-point correlation functions between the left-handed currents and the 
weak operators, averaged over the spatial volume, now read [15]: 

Tv[T a T b ]C(x ) = Jd 3 x(jg(x)J b (0)) , (2.16) 

[Cn] a \x ,y Q ) = Jd 3 xjd 3 y(jg(x)O R (0)J b (y)), (2.17) 

where the index R refers to the representation. 

On the xPT side, the operator corresponding to Eq. (2.15) becomes, at leading order in 
the momentum expansion, 

(2.18) 

The two-point correlation function C(xq) is defined (apart from contact terms) by 

Tr [T a T b ] C(x ) = j d 3 x ( J£(x)J b (0)) , (2.19) 

and the three-point correlation function we are interested in, reads (again apart from contact 
terms) 

[CR] a \x ,y ) = J^J d 3 y (j a (x)0 R (0)Jo(y)) ■ (2-20) 

Our task is to compute the objects in Eqs. (2.19), (2.20) under certain circumstances, to be 
specified in the next Section. 

3. Regimes of chiral perturbation theory 

Given a fixed spatial extent L » 1/F of the box, several different kinematical regimes can 
be identified in xPT, leading to various computational procedures [45]. The situation is 
summarised in Fig. 1. We will here be interested in the p- and e-regimes; the 5-regime 
(corresponding to small but elongated boxes) is also relevant in principle, but quite tedious 
to handle in practice [45], and thus preferably avoided. 

3 The convention in Eq. (2.15) differs by a factor i from that in Ref. [19], but agrees with the convention of 
Refs. [17, 20]. We use this "unphysical" convention since it removes a number of unnecessary overall minus 
signs from the xPT predictions. 
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J a = Tr 
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Figure 1: The different regimes of chiral perturbation theory, given a fixed spatial extent L of the 
box, according to Ref. [45]. Here T is the temporal extent of the box and M the pseudoscalar mass. 
It is assumed that L> l/F. 

3.1. p-regime 

In the p-regime, the quark mass is large enough to ensure that 



It follows from this condition that the Goldstone field £, defined through U = exp(2i£/F), 
behaves effectively as a small quantity, and can be expanded in. Chiral corrections are 
obtained as an expansion in {M/Ff and l/(FL) 2 , where M 2 = 2mE/F 2 . The power- 
counting rules in this regime count both of these expansion parameters at the same order: 



where p is assumed small, p <C F. The temporal extent T can in principle be small or large, 
as long as T > 1/p. Of course, it is also possible to send L — > oo in the p-regime expressions. 
The situation is illustrated in Fig. 1. 

Inserting the Taylor-series of U into Eq. (2.9), the propagator becomes 



mT,V > 1 . 



(3.1) 



M ~ p , L ~ — 



(3.2) 



P 




(3.3) 



where 




(3.4) 
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and V = TL 3 is the volume. Here we have also set = 0, as is usually done in the p- 
regime. In the unquenched case, E(x; M 2 ) = G(x; M 2 )/N{, but we keep everywhere E(x; M 2 ) 
completely general. The reason is that then the form of Eq. (3.3) is general enough to contain 
also the propagator of the replica formulation of quenched chiral perturbation theory [46, 13]. 

For future reference and as an example of a NLO result in the p-regime, we consider the 
two-point correlation function in Eq. (2.19). The result can be written as 



C(x ) = 



+ 



1 



E(0;M 2 



W(0;M») + 5 £ ■ 



F 2 

where (for \xq\ <T) 



F 2 

8M 2 



F 2 



M 2 P(x ) - 



Nf dG(0; M 2 
_ 



(7V f L 4 + L 5 - 2iV f L 6 - 2L 8 ) 



d 



dM 2 



dT 



M 2 P(x ) 



P(x ) = |d 3 xG(x;M 2 ) = i^ 



e ip x _ cosh[M(T/2 - |x |)] 
2Msinh[MT/2] 



T ^ vl + M 2 

PO u 



while 



G(0; M 2 ) = Goo(M 2 ) + G V {M 2 ) , 
where Goo(M 2 ) is the infinite-volume value, 4 

d d p 1 
(2vr) d p 2 + M 2 ' 



Goo(M 2 ) = J 



and the (finite) function Gy(M 2 ) incorporates all the volume dependence [47], 1 



G V (M 2 



(47T) 



dA 



-AM 2 



E( 1 -«) exp [-4T( r2n ° +L2|n|2 ) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



For MV i ^> 1, the finite-volume effects are exponentially small, and we can set Gy = 0. 



3.2. e-regime 

In the e-regime, the natural dimensionless variable is fx = mT,V. The power counting rules 
are now 

mS~e 4 , L ~ - , T ~ - , (3.10) 

where e is assumed small, e <C F. Of course, it is also possible to send m — > in the e-regime 
expressions. Hence the parameter /i is parametrically of up to order unity. In this regime, 
the Goldstone boson zero-mode Uq, defined by writing U = exp(2i£/F)Uo, where £ has non- 
zero momenta only, dominates the dynamics, and needs to be treated non-perturbatively. 

4 The divergence of Goo{M 2 ) for d « 4 cancels against those in the Li's [43], cf. Eqs. (B.17), (B.31). 
5 In Ref. [47] the function Gv was denoted by gi. 



7 



Consequently, gauge field topology plays an important role [48], and it is useful to give the 
predictions in sectors of a fixed topological charge v. 

As an example, the two-point correlation function C(xq) of Eq. (2.19) becomes [49, 14, 19] 



C(x ) 



2T 



,,N t ( fa T 2 k 00 \ , 2TV 



(3.11) 



where xo = xq/T, and the constants j3± and &00 are related to the (dimensionally regularised) 
value of 



ip-x 



P 



2 ' 



by 



G(0) = 



V ' dT v 7 y 



Introducing p = T/L and 



(3.12) 



(3.13) 



(3.14) 



where S"(x) is an elliptic theta-function, <S(x) = J2^=-oo ex P( — itxn 2 ) = #3(0, exp(— 7rx)), a 
numerical evaluation of these coefficients is possible through (see, e.g., Refs. [47, 49]) 



n^O 



sinh (7rp|n|) 



(3.15) 
(3.16) 



Furthermore, <T v (fi) = Af~ d{lndet[J^+j_j(/x)]}/d/Lt, where the determinant is taken over an 
iVf x iVf matrix, whose matrix element is the modified Bessel function I v+ j^i [50, 48]. 
The function h\{r) appearing in Eq. (3.11) reads (for |r| < 1) 



Mr) ^ - 2 



1 

T| ~ 2 



1 

12 



(3.17) 



3.3. Further remarks 



In the following, we carry out computations according to the p and e-countings as outlined 
above. Other recent work for related observables has made use of the p-regime, with T 3> 
L [51, 52, 53]. There have also been extensive NLO computations at infinite volume [54], 
which is a special limit of the p-regime. 

Note that if 1/FL <C 1 as our power-counting rules assume, and we consider an observable 
that is independent of the topological charge v, then the e and p-regimes should in principle 
be continuously connected to each other (cf. Fig. 1). Concretely, for ML <C 1 and T ~ L, 
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Figure 2: The NLO graphs for C27 in the p-regime. Lines denote meson propagators, an open 
square the left-handed current, an open circle the weak operator, and four-point interactions with no 
symbol and with a closed circle the "kinetic" and "mass" terms in the chiral Lagrangian, respectively. 
Diamonds indicate QCD and weak interaction 0(p ) low-energy constants. 

Eq. (3.5) goes over into Eq. (3.11) with /x 3> 1, in which limit the dependence of Eq. (3.11) 
on v disappears. Whether such a crossover takes place in practice remains to be inspected 
for each observable separately, and gives some feeling concerning the convergence of the xPT 
computation, i.e., whether 1/FL <C 1 is satisfied. 

4. The AI = 3/2 operator 

We now address the determination of g27, considered previously in the e-regime [19]. 



4.1. p-regime 

The graphs entering the computation of Eq. (2.20) at next-to-leading relative order in the 
p-regime are shown in Fig. 2, with the weak operator O27 to be taken from Eq. (2.11). The 
result can be written in the form 

[C 27 ] ab (^o, yo) = A$ [c(x )C(y ) + V 27 (x , y )] , (4.1) 

where (for N v = 3) 

o o 

As an example, choosing kaon and pion type currents, we could take 

T^ = 6 iu 5 js & J$ = u 7o P-s, (4.3) 
T^ = 5 ld 5 ju & 4 = d lQ P-u, (4.4) 
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Figure 3: The function [C 27 }^ orm {-T/3,y ). The parameters are: N f = 3, F = 93 MeV, L = 2 fni, 
T/L = 2, A = 1000 MeV. 



and then 



a ab 
^27 



(4.5) 



Given the result for C(xq) in Eq. (3.5), the only further missing ingredient in Eq. (4.1) is 
T^27(xo,yo)- We obtain 



m , F 2 (M 2 d 2 G(0;M 2 

V2 7 (x ,y ) = -—\— dM 2 dT 



+ M 



2 dG(0;M 



2^ r 



dT 



P(x )+P(y ) 



+ 



+2M 4 G(0; M 2 )P(x )P(y ) - -M 4 P(x - y ) B(x ) + B(y ) 



+ 



+M 4 £dr P'(t - x )P'(t - y )B(r) j , 
where the new object B(xq) is defined as (for |xq| <T) 



B(x ) = J d 3 x G(x;M 2 



L 3 S 



cosh[£(T/2 - |x | 
2£sinh[£T/2] 



(4.6) 



(4.7) 



E=^M 2 +p 2 



The expression in Eq. (4.6) is, as such, ultraviolet divergent: in dimensional regularization 
in d = 4 — 2e dimensions, the third and the last terms on the right-hand side contain poles 
in e. Denoting A = — l/327r 2 e, we can write 



V 27 {x ,y Q ) = V r 27 (x ,y ) + F 2 \[^M 4 P'(x )P'(yo) - M 6 P(x )P(y ) 



(4.8) 



where T>2j(xo,yo) is finite. The divergences get cancelled against the 0(p 4 ) low-energy con- 
stants related to weak interactions, as shown in Appendix B. As there are a large number 
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Figure 4: The values of [C 27 ]£orm(- r A T /^)- Th e parameters are: N f = 3, F = 93 MeV, L = 2 fm 
(left), L = 4 fm (right), T/L = 2, A = (500 - 2000) MeV. 



of them, however, it is sufficient for our purposes here to note that the 0(p 4 ) low-energy 
constants amount to cancelling the 1/e-divergences in the result and replacing the corre- 
sponding MS scheme scale parameter Ji by two different physical scales, A for the coefficient 
of P(x )P(y ) and A' for the coefficient of P'(x )P'(yo)- 

For practical applications, it is convenient to normalise the three-point correlator by divid- 
ing with two two-point correlators: 



The function 1Z27(xo,yo) is then trivially obtained from Eqs. (4.6) and (3.5); in Eq. (3.5), it 
is even enough to keep the leading order contribution only, since V 2 7 (xq, j/q) gets generated 
only at NLO. 

As an example, the function [C27]no rm (— T/3, yo) ^ s plotted in Fig. 3 as a function of yo, for 
the index choice in Eq. (4.5) (solid line). The values of [C 27 ]f orm (-T/3,T/3) are shown in 
Fig. 4, as a function of ML (the region bounded by solid lines). In these plots, the effects of 
the weak LECs have been collected to a single scale A = A' appearing inside the logarithms, 
and the scale has been varied in a wide range, to indicate the size of the uncertainty related 
to the unknown higher order LECs. 

We would like to stress at this point that the p-regime results are parametrically valid only 
in the range ML>1/FL: for generic observables, the contributions of the Goldstone zero- 
modes become dominant if this inequality is not satisfied, and need to be resummed, leading 
to the rules of the e-regime. It turns out [20], however, that in the normalised observable 
[C27] norm ( x 0) 2/o ) that we have considered here, the contributions from the Goldstone zero- 



[C27] ab (x ,yo) _ Aa6 r i , V 27 {xQ,yQ 
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modes cancel out at this order. Therefore the result can in fact formally be expanded as a 
Taylor-series in (ML) 2 , with the zeroth order term agreeing with the result of the e-regime 
(see below). Still, one has to keep in mind that the Taylor-expanded result only needs to 
reproduce the correct mass dependence in the range ML > 1/FL. 

Let us finally briefly touch the conventional limit of large volumes. We assume xq = 
— \xo\,yo = 1 2/0 1 ^ such that the charges are on opposite sides of the operator. Then P(xo) = 
exp(-M|x |)/2M and P'(x )P'(yo) = -M 2 P(x )P(y ). In other words, the distinction 
disappears between the two structures getting contributions from the higher order LECs (cf. 
Eq. (4.8)), just as would happen if a partial integration could be carried out with respect to 
the position of the weak operator. Consequently, only a single combination of LECs appears, 
and the corresponding effects can be collected into a single scale A. We obtain 



^■27(^0,yo) 



M 2 



where 



(4vrF) 2 



(x 



3 In + 2 - e - 2M ^l0(2M|xo|) - e~ 2M ^\ ^M^) 



s f 

Jo 



dz z^e 



l + z 



+ 



2 + z l + z 



- 2 



(4.10) 



(4.11) 



The xo and yo-dependences in Eq. (4.10) are very small in practice. As seen in Fig. 4 (dotted 
line) , one needs to go to volumes as large as ML > 5 in order for the simple infinite- volume 
approximation to be accurate for this observable. 6 



4.2. 



e-regime 



The e-regime results for T>2t{xq, y$) were derived in Ref. [19] but, for completeness and future 
reference, we briefly reinstate them here. For V27 in Eq. (4.1) one obtains 



^27(^0,2/0) 



2T 2 



1 r^JG(o), 



(4.12) 



and, using Eq. (3.13) as well as the leading-order part of Eq. (3.11), the ratio in Eq. (4.9) 
becomes 



^27(^0, yo) 



{FLf 



p 2 Pi- pk 00 



(4.13) 



where p = T/L. Note that this result is independent of the topological charge v, although 
computed in a fixed topological sector. 



The e-regime prediction for the function [C27]^ orm (— T/3, yo) is plotted in Fig. 3 as a function 
of yo, for the index choice in Eq. (4.5) (dashed line). The values of [C27]norm( — T/3, T/3) are 
shown in Fig. 4, as a function of p (dashed line). 

6 Note that finite-volume corrections depend on the observable in question; in particular, the finite-volume 
effects that we find are much larger than those in typical two-point correlation functions. 
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4.3. Further remarks 



Let us inspect Fig. 4(left), around the region ML ~ 1.5, or fj, ~ 2.0. Moving to smaller values 
of fi, the e-regime becomes more accurate, while at larger ML, the p-regime should be the 
correct procedure. But which result represents better the truth at this intermediate point, 
where both countings are in principle parametrically applicable? 

Let us note that for the semi-realistic parameters used in Fig. 4(left), 1/FL ~ 1.1. There- 
fore, the parametric rules we have assumed are at best satisfied by a narrow margin. Con- 
sequently, higher order corrections in xPT can be important. In the absence of an explicit 
computation thereof, it remains to be inspected phenomenologically which of the predictions 
reproduces better the volume and mass dependences of the simulation results in this regime. 

We end with a small remark on quenching. Employing the replica formulation [46, 13] of 
quenched chiral perturbation theory [55, 56], the only changes with respect to the unquenched 
situation are that we need to replace the propagator of Eq. (3.3) through 

£(*;M 2 ) = ^G(x-M*) + ml -* M2 H( X -M*) , (4.14) 
1 x-^ e ip ' x 

H( X ;M>) = ( 2 + M2)2 , (4.15) 

where new parameters related to axial singlet field, mQ/2N c , a/2N c , have been introduced; 
and take Nf — > at the end of the computation. Given that our results for [C27]^ rm (xo, yo) 
are completely independent of iVf and of the function E(x; M 2 ), however, there is no change 
with respect to the unquenched theory for this observable [19]. 



5. The AI = 1/2 operators 

In the case of the AI = 1/2 transitions, two operators with the right symmetries appear in 
Eq. (2.10). This means that if we have measured some correlation function on the QCD side, 
with an operator Og transforming in the octet representation, then this is to be matched to 
a linear combination of correlation functions on the side: 

Jd\ ld%^JS(x)h 8 Os(0)4(y))^gs[Cs] ab (xo,yo)+g's[Cs} ab ^o,yo) , (5.1) 

where /ig is the Wilson coefficient, and gs ,g' s are the partial contributions from /ig Og to the 
corresponding LECs. We thus have to consider two different classes of correlators on the 
xPT side, in order to be able to disentangle the coefficients of these operators. 

5.1. p-regime 

For the operator Og of Eq. (2.13), the graphs entering the computation of the correlation 
function in Eq. (2.20) are the same as in Fig. 2, and the correlation function has the same 
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form as in Eq. (4.1): 



[Cs] ab (x ,y ) = Af C(x ) C(y ) + V 8 (x , yo) 



, ab 



where 



Af = -{T a ,T b } ds , 



(5.2) 



(5.3) 



and the function C{xq) is still given by Eq. (3.5). For the matrices T a ,T b in Eqs. (4.3), (4.4), 
the group theory factor evaluates to 

Af = \ . (5.4) 

The function 2?g(xo,2/o) reads 



^8(^o, yo) = -^^27(^0,2/0) + 



+ F 2 M 2 ^±1 



G(0;M 2 )-2E(0;M 2 ) P'(x )P'(y ) + 



1 



1 



+ 



+M z P{x - yo) [B(x ) + B(y ) - -B{x ) - -B{y 
+P'(x - y ) [&{yo) + Bo(yo) + \b'{x q ) - B'(x ) - £ (x ) - ^B'(y ] 



+M 2 [ dr 
Jo 

The new objects appearing here are defined as 

B(x ) = fd 3 xG(x;M 2 )E(x;M 2 ) , 



M 2 B(t) - 2M 2 B(t) - B 00 (t) P(t - x )P(t - y ) 



Bo 
Boo 



d G(x; M 2 )E(x; M 2 ) - G(x; M 2 )d £(x; M 2 ) 



(so) = y«i 

(x ) = y"d 3 x d 2 G(x;M 2 )E(x;M 2 ) - G(x; M 2 )d^E(x; M 2 



+ 
(5.5) 

(5.6) 
(5.7) 
(5.8) 



We recall that in the unquenched theory, E(x;M 2 ) = G(x; M 2 )/N{, and B(xo) thus agrees 
with B(xq)/N{ as defined through Eq. (4.7), while Bq(xq), Bqq(xo) vanish. 
Eq. (5.5) again contains divergences: in the unquenched theory, 



V 8 (x ,yo) = Vl(x ,yo)+F 2 \[(--^±_)M i P>(xo)P'(yo) + 

1 )M 6 P(x )P(y ) 



N v - 2 iV v + 2 - 
4 + 2Nf J' 



(5.9) 



where A = — l/327r 2 e, and Dg(xo,2/o) is finite. The cancellation of these divergences against 
the 0(p 4 ) LECs is demonstrated in Appendix B. 
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Figure 5: Left: the function [C 8 ]^ rm (-T/3, y ). Right: the function [C' 8 ]^ orm {-T/3, y ). The pa- 
rameters are: N { = 3, F = 93 MeV, L = 2 fm, T/L = 2, A = 1000 MeV. 



Following Eq. (4.9), it is convenient to define a normalised correlation function by dividing 
with two current-current correlators, and we thus obtain 



[C 8 ] afc (x , y ) 
J,lonu ' ~ C(x )C(y ) 



[C 8 ]l m (* ,2/o) = 



A 



ab 



1 +K 8 (xo,yo) 



(5.10) 



Again, it is enough to use the leading order forms for the functions C(xq), C(yo) in the 
definition of TZs(xo,yo), since Dg(a;o,yo) gets generated only at NLO. 

In the infinite- volume limit, the distinction between the various types of divergences in 
Eq. (5.9) disappears, as before. Collecting the corresponding LECs to a single scale A, we 
obtain (in the unquenched case) 



Ks{xo, 2/o ) = -^^-27(^0,2/0) + — 2 



2 \ M 2 



A 2 

2 In — + 1 + 



Nf J (4ttF) 2 [_ M 2 
+e~ 2M ^~(2M\x \) + e' 2M ^E(2M\y 



where 7?.27(^0; 2/0 ) is from Eq. (4.10), and 



H(x) = 



dzz2 e 2:2 

1 + z 



2 + z 1 + z 



+ 2 + 4z 



, (5.11) 



(5.12) 



For the correlator C' 8 the graphs are the same as in Fig. 2 except that, for a vacuum angle 
= 0, the weak operator 0' s only couples to an even number of Goldstone modes. The result 
is now of the form 

(5.13) 



[C' 8 } ab (x ,yo)^A a 8 b V' 8 (x ,y ), 



ab -7-)/ , 
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Figure 6: The function [C 8 ]^ rm (-T/3,T/3). The parameters are: JV f = 3, F = 93 MeV, L = 2 fm 
(left), L = 4 fm (right), T/L = 2, A = (500 - 2000) MeV. 

where 

£>s(>o,yo) = 



F 4 j 


1 - 


"2" I 


NfM 4 


' 2F2 J 


2M 


-( 


F 2 


■/ 



jV f G(0;M 2 ) | E(0;M 2 ) m2 d 



M 2 P'(x )P'(y ) 



+ 



F 2 F 2 dM 2 

dr [>( r _ x )P'(r - y ) + M 2 P(r - x )P(t - y )\ B{t) 

N f M 2 dG(0; M 2 ) 



dTP / {T-xo)P'{T-yo)B(T) 



(5.14) 



2F 2 T dM 2 dT 
Separating the divergent parts, we get (in the unquenched case) 

V 8 (x ,yo) = V'£(x ,y ) + F 2 \[(-^ + ^)M 4 P\x )P'(yo)-^M (i P(x )P(y ) + 

(5.15) 



+ ^ M6 d^( p/ ^ p/ ^ 

where 2?g r (x ,?/o) is finite. The cancellation of divergences is demonstrated in Appendix B. 

If we want to disentangle the dependences following from the operators Og and 0' s in 
a given lattice measurement, we are lead to compare the contributions from 0' 8 with the 
normalised correlation function in Eq. (5.10). Therefore, we define 



[CsUJ^yo)= c{xo)c{yo) 



(5.16) 



Treating UV-divergences and higher order LECs as before, the correlation functions 



[ c s]norm(- T A 2/0 ) and [Cg]no rm (-r/3,y ) are plotted in Fig. 5 as a function of y (solid 



1 ab 
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Figure 7: The function [C' 8 ]l b olm (-T/3, T/3). The parameters are: N f = 3, F = 93 MeV, L = 2 fm 
(left), L = 4 fm (right), T/i = 2, A = (500 - 2000) MeV. 



lines). The two correlators are observed to have a rather different dependence on yo, so it is 
in principle possible to disentangle their contributions in a given lattice measurement. The 
values of [C 8 ]^ rm (-T/3, T/3) and [C^ rm (-T/3, T/3) as a function of ML are illustrated in 
Figs. 6, 7, respectively (regions bounded by solid lines). 

It is important to stress that, for ML — > 0, the correction of relative order 1/F 2 in Eq. (5.14) 
diverges as ~ l/i^M 2 ^. This indicates in a concrete way that the p-regime computation is 
no longer reliable for ML <C 1/FL, and we need to turn to the e-regime. 

Let us again end by commenting on the conventional limit of large volumes. Assuming 
xo = —\xo\,yo = 1 2/0 1 5 and inserting the unquenched value of E(x;M 2 ), we obtain for the 
normalised case 



-{T a ,T b } ds \l + 



M 2 



(4ttF) 2 



-N t [ 1 + ln 



A^ 

M 2 



2 , A 2 
— In — ^ 

N { M 2 



where 



A(x) 
T(x) 



_^Vf | e _2J»#|xo| A(2M|x |) + e- 2M l^lA(2M|y |)} 
+-^{e~ 2M ^r(2M\x \) + e- 2M l y °lT(2M|y |)} 

f 

Jo 

f 

Jo 



+ 



dz z^e 



dz z2e~ 



: V2Tz~ 
l + z 

: V2Tz 
l + z 



2 + z 
1 



2 + z + 1 + z 



, (5.17) 

(5.18) 
(5.19) 
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Figure 8: The leading-order graphs for [C' s } ab (x ,yo) in the e-regime. An open square denotes the 
left-handed current, an open circle the weak operator, and a filled circle a mass insertion. 

This result is plotted in Fig. 7 with dotted lines. We observe again how only values ML > 5.0 
guarantee that finite-volume effects are small for our three-point observables. 

5.2. e-regime 

We finally move to the e-regime. For Cg the graphs are the same as for C27, as depicted in 
Fig. 3 of Ref. [19]. The correlator retains the form in Eq. (5.2), with A 8 b from Eq. (5.3), C{xq) 
from Eq. (3.11), and T>27(xo,yo), appearing as in Eq. (5.5), from Eq. (4.12). The order of 
magnitude of the leading term in [Cg] afe (xo, yo) is 0(e 2 ), and the NLO term is 0(e 4 ), while the 
terms beyond T>27(xq, yo) in Eq. (5.5) are formally 0(e 6 ), so that the corresponding graph 
(the sixth in Fig. 2) can be ignored in the e-regime. Therefore, all information is in the 
e-regime version of T>27(xo, Vo)- To be explicit, the normalised form of Eq. (5.10) becomes 



ML>o,yo) = Af 



(FL) 



(5.20) 



Let us stress, in particular, that [Cg]^ rm (xo, yo) is independent of topology and quenching 
at this order, just like [C27]^ orm (xo, yo)- (However, as discussed at the end of Appendix C, 
quenching does lead to the appearance of additional LECs [32] that need to be disentangled.) 

Let us then address C' 8 . Given that the e-regime computation is to be carried out at fixed 
topology, the operator 0' 8 needs now to be considered in the full generality of Eq. (2.14), i.e., 
with a non- vanishing vacuum angle 6, unlike in the p-regime. Therefore 0' 8 can also couple 
to an odd number of Goldstone fields. On the other hand, it is easy to see that the tree-level 
graphs (cf. Fig. 8) are already of order 0(e 4 ). Comparing with Cg, it is therefore enough to 
restrict to the leading order. We find 

[C'sfi^yo) = ^{{T a .T b } ds a u ^)h' 1 (xo)h' 1 (y ) + 

+ [T a ,T b ] ds ^{ti 1 (x - yoMi&o) + h[(y )] + hx(y ) - /*(£,,)} j , (5.21) 



where hi is from Eq. (3.17). The corresponding normalised form reads 

4T^ 
-pi 



ab ¥T 2 , ah 

[Cslnorm 0*0,2/0) = [4] (^0, ^o) ■ (5.22) 
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Figure 9: The graphs contributing to [/Cg] a (a;o). The notation is as in Fig. 8, with additionally a 
cross denoting a "measure term" (cf. Ref. [9]). 

This result is to be used in combination with Eq. (5.20), in order to disentangle the two terms 
on the right-hand side of Eq. (5.1). 

Unlike Eq. (5.20), the expressions in Eqs. (5.21), (5.22) get modified in the quenched theory, 
because they contain Goldstone zero-mode integrals. Proceeding as in Ref. [16], we find 



where I u , K v are modified Bessel functions. Note that Eq. (5.23) could also be obtained from 
Eq. (5.21) by just naively setting N{ — ► and replacing a v — > a qu . 

Since the functions [Cg] ab (xo, yo) an d [C' 8 ] a (xo,yo) are not identical, a precise measurement 
of the time-dependence of the correlation function of the left-hand side of Eq. (5.1) would in 
principle make it possible to disentangle the contributions to gs,g's- in particular, as shown 
in Fig. 5, any dependence of the correlation functions on v arises at this order through the 
operator 8 . In practice, however, the problem emerges that it may not be easy to obtain 
such a high accuracy that the two LECs could reliably be determined from a single observable. 
Therefore, it may be beneficial to define another probe as well, such that the LECs can be 
disentangled with better confidence. We now show how this can be done. 

5.3. Direct determination of g' 8 

In order to determine g' 8 , we consider the correlator 



where the subscript q refers to the quenched theory, and [57] 

a qu (n) = n\l v {ii)K v {n) + I v+x (ji)K v -.i.(n) + ~ , (5-24) 





(5.25) 



on the side of QCD, and correspondingly 




(5.26) 
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on the xPT side. Note that this correlation function is not available in the conventional 
p-regime setup (i.e. with 6 = 0), because it is odd in charge conjugation. 

We have computed both [/Cs] a (xo) and [/Cg] a (x ) at NLO in the e-regime. Parametrically, 
the orders of magnitude of the LO and the NLO graphs are 0(e 4 ) and 0(e 6 ), respectively. 
We find, however, that at this order [/Cs] a (xo) vanishes exactly, like in the p-regime. 

On the other hand, [IC' 8 } a (xo) does not vanish. The graphs are shown in Fig. 9. We find 



vF 2 



ds ' 



V 



h'i(x ) 



1 + 



1 

iVf 



G(0) 
F 2 



+ h' 2 (x ) 



2T 2 



1 

iVfJ 



where G(x) is from Eq. (3.12), and (for |r| < 1) 



r 2 (M 



I) 2 - ± 
1 30 



The result is illustrated in Fig. 10, after normalisation through 

L 3 [/C^(* ) 



C(x ) 



= ^Xl^-slnorm^o) 



Repeating the same steps in the quenched theory, we find 

vF 2 



KQxo) = -T a ds - 



V 



h'i{xo) 
2T 2 

'fW 



1 + 



a G(0) ml H(0) 



2N r F 2 2N r F 2 



+ 



(^M + ^r) h> 2^o) - ^h' 3 (x ) I , (5.30) 



(5.27) 



(5.28) 



(5.29) 
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where H(x) and H^{t) (for |r| < 1) are defined through 



(p2)2 



^ 3(T) 55 720 
The value of #(0) is given by [47] 



t 2 (\t\-1) 2 (r 2 -\r\-^)+i- 



42 



H(0) 



02 + 



-2e 



(4vr) 2 [e 



+ ln(/iV 1 / 2 ) +l + 0(e) 



where In// = ln/r + ln47r — 7_e, and (with a from Eq. (3.14)) 



ft = 



(4vr) 2 



"o(pSP 4 )+a_ 2 (p 4 ,P 4 ) - 2 ~ m ( 47 + 7s 



(5.31) 
(5.32) 



(5.33) 



(5.34) 



The UV-divergence in Eq. (5.33) is cancelled by £ (cf. Eq. (2.14)), which is to be treated as 
a bare parameter in the quenched theory [58] . 

To summarise, we now have a method to disentangle the two contributions related to the 
LECs gs,g' 8 : by considering [-^slnormO^o)) we can first match for g' 8 . Then the corresponding 
term can be subtracted from the right-hand side of Eq. (5.1), and we are able to determine gs- 
As illustrated in Fig. 5, a cross-check is that the dependence on v should have disappeared. 



5.4. Further remarks 

The remarks that can be made on the convergence of the e and p-regime computations of Cs 
and C' 8 are largely the same as for C27 in Sec. 4.3. Indeed, for 1/FL <^ 1, there could be a 
non-vanishing overlap, i.e. a regime where both the p-regime and the e-regime expressions 
are valid. For the more realistic case 1/FL ~ 1, on the other hand, this is unlikely to happen. 
It would be tempting to read from Fig. 7 that the p-regime expression works in the range 
ML > 2.0, and the e-regime expression in the range fi < 2.0, but whether this is really the 
case remains to be seen once a comparison with lattice simulation results is available. 

Concerning quenching, let us stress that the correlation function [Cs]^ rm (xo, yo) is de- 
termined in the e-regime by the same function 7^27(^0 > yo) as [C27]norm( x o> 2/0 )> and is thus 
insensitive to quenching at the present order. At the same time, the correlation functions 
[ C s]norm( x o,yo) and [/C 8 ]^ rm (x ) do get modified. 

An important point however is the different relevance of the quenched ambiguities of 
Ref. [32] in the two regimes. 7 In general the quenched theory contains spurious operators 
with new LECs. Some of these originate from the fact that Nf ^ N v , a case that is considered 

7 We refer here to the ambiguities at the level of the chiral Lagrangian. We assume always that the weak 
effective Hamiltonian at the quark level contains an active charm so that no "unphysical" operators appear 
in the Operator Product Expansion at the order in the Fermi constant at which we are working. 
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in detail in Appendix C, while others are related to the couplings of the axial singlet field 
that cannot be integrated out in the quenched limit. The latter modify the terms that in 
the full theory would be divergent in the limit Nf — > 0. We indeed confirm a rather messy 
situation in the p-regime, where many new couplings enter; thus we have not carried out a 
systematic study of all quenching effects in our observables in this regime. On the other hand, 
the quenching ambiguities are reduced to a minimum at the NLO in the e-regime, with appar- 
ently only one spurious octet LEC contributing to [CslnormC^Oi Vo)- Therefore certain octet 
couplings can be determined by matching the lattice simulation results to [Cs]^ orm (xQ,yQ). 
We elaborate on this issue in more detail in Appendix C, particularly around Eq. (C.23). 

6. Conclusions 

We have addressed in this paper the determination of the 0(p 2 ) LECs of the chiral weak 
Hamiltonian. As probes we have used the three-point correlation functions between the weak 
operators and left-handed flavour currents. We have computed the three-point correlation 
functions up to next-to-leading order in chiral perturbation theory, both in the e and in the 
p-regimes, for all three operators that appear in the SU(3) chiral weak Hamiltonian. 

While the determination of the LEC 527 , which fixes the AI = 3/2 amplitude of the 
weak decays K — > tttt as well as the kaon mixing parameter Bk in the chiral limit, appears 
straightforward, the determination of the LECs fixing the AI = 1/2 amplitudes is more 
demanding in several respects. Even restricting to the idealised case of full QCD at large 
volumes, there are two operators with the same flavour symmetry, while only the coefficient 
of one of them, g$, contributes to the physical kaon decays [3, 4]. Therefore it is important 
to come up with a setup which makes it possible to remove the contamination from the other 
operator in a lattice measurement of the type that we have considered. 

We have shown here that this challenge can be met by going to the e-regime. The two 
operators contribute in very different ways to a given three-point correlation function, one 
leading to a topology-dependent and the other to a topology-independent result. Moreover, 
we have found a two-point correlator that is only sensitive to the "unphysical" LEC and 
can be used to fix it. Therefore, it seems possible in principle to disentangle the physical 
coefficient gs from lattice measurements in the e-regime. 

By comparing the e-regime results with p-regime results in a finite volume, we have also 
speculated on the regimes of validity of the two approaches. It appears that for semi-realistic 
lattices with a spatial extent of about 2 fm, the e-regime approach might be applicable for 
H = mT,V < 2.0 and the p-regime for ML > 2.0. In any case, the conventional infinite-volume 
formulae are accurate (with errors below 10 - 20%) only at ML > 5.0. 

Finally, we have briefly addressed the effect of quenching on the determination of the 
AI = 1/2 observables. New unphysical couplings are in general expected in the effective 
chiral theory with respect to the unquenched situation. In the e-regime we find, however, that 
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the contamination from these new couplings is minimal at NLO: only one additional coupling 
enters our predictions, and we have shown that it is in principle possible to determine the 
quenched g$ in spite of these quenching artifacts. 
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Appendix A. Irreducible representations of the valence group 



For completeness, we reiterate in this Appendix the main formulae related to irreducible 
representations of the valence group SU(iV v ), relevant for the operators appearing in the 
weak Hamiltonian. We follow the tensor method discussed, e.g., in Ref. [59]. 

Like in the main body of the text, we make a distinction between the valence group 
SU(iV v ), used to classify the weak operators, and the full flavour symmetry SU(iVf). The 
indices f,s,u,v,f,s,u,v,f,s,u,v are assumed to take values in the valence subgroup only. 
We denote by Ofsm a generic operator transforming under N* <g) N* <g> N v (g) N v of SU(iV v ), 
and by <9~ s one transforming under N* <S> N v . 

We define the projection operators 

(-Pi )fsuv;fsuv = ~^(.dffdss ~\~ &5fs$sf){$uu$vv ~\~ <J$uv$vu) , (A-l) 
(^°2 ^rsuv;rsuv = ^ff^ss^uu^vv ~\~ TTT ; « \/»r I \{^ru^sv ~\~ &8rv8su)8f u 5sv 

(N v + 2a)(N v + a) 
( $fu $ss $vv 8f-u + 5 sv 5 rr 8 uu S S v ~\~ o'S rv 5 ss 5 U q)5f U + (j5 su 8 rr 8 vu 5 sv ) , (A. 2) 



N v + 2a 



{P?i)ru;fu — 8 r f5 uu 8 ru 8 ru . (A. 3) 

In addition, P v is defined to project from SU(iVf) to SU(7V V ). The operators denoted by O27, 
08 and 0' 8 can now be defined as 



[027}fsuv = ( P 2 P l)fsuv;fsuv Ofsui) , (A.4) 
[^8 \fu ~ (^3) fu;ru(Rl )fsus;fsuv Ofsuv J (A-5) 
[O'sU = (P3)ru;fuO'ru • (A.6) 

Note that the contraction over I in Eq. (A. 5) goes over valence flavours only, and that 
additional octet operators (0%) can appear already at the leading order when N v 7^ Nf. 

Instead of a generic operator Ofsm-, practical computations of the type in Ref. [60] involve 
certain factorised forms, like 

[OiU- fflS = (Q)uf(R)vs, (A.7) 
= (Qhs(R)w , (A.8) 
[0 3 ]_„ s = 5 U r(R) v - s , (A.9) 

PaVsuv = Sus(R)vf, (A.10) 
[O^rsu-v = SurSvs, (A.ll) 
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Then projections of the types in Eqs. (A. 4), (A. 5) produce 



[PgPfOi]—-- 


^rsuv (.Qi P) J 




— PfvkQi p) > 




= aS^(Q,R), 




= aT f %(Q,R) 


[^2 ^1 O&rsuv 


= 0, 


Ip 3 p?o 3 U 


= U~ a (R) , 




= 0, 


[P 3 Pi a o A U 


= <tU^(R), 


[P2 Pi ^5 ] fsuv 


= 0, 


[PsPi a 5 U 


= 0, 


[PiPlOelrsuv 


= 0, 


\P,PiO,] m 


= 0, 



(A.13) 



where (introducing the notation Tr v (...) = Tr (P v ...)) 

S?suv(QiP) = 



QufPvs PufQvs ^(.QusPvf PusQvr) 

[ (QP V R + RP v Qhs + a(Q us Tv v (R) + i^ s -Tr v (Q)) 
{QP V R + RP V Q) W + <t(q w Tt v (R) + R w Tt v (Q)) 
(QP V R + RP v Q)v- s + a^ s Tr v (R) + R^Tr y (Q)) 
(QP V R + RP v Q)uf + cr(Q af Ti v (R) + R af Ti v (Q)) 
Tr v (QP V R) + aTr v (Q)Tr v (R) 



N Y + 2a 
N v + 2a 
N v + 2a 
N v + 2a 



(N v + a)(N v + 2a) 



T r UQ,R) 



a{QP v R + RP v Q)uf + Q Sf Tr v (R) + R Q , Tr v (Q) 



2iV v 



U? U (R) = -(N v + 2a) 



Tr v (Q)Tr v (R) + aTr v (QP v i?) 

5, 



R uf - ^ Tr v (^) 



+ 

(A.14) 

(A.15) 
(A.16) 



Considering, in particular, the operators 

(A.17) 
(A.18) 

(A.19) 
(A.20) 

which appear in the computations of Fig. 2, and choosing the indices that appear in the 



^■rsuv — Pus^vr (j^us{P ->P }vr ^vr{P ->P }usj j 

aI-L — T^- T-1 + lT a TH I- r) lT a T b T-- 



TSMO 

A<1L 

rsnu 
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physical operators O27 and 8 , we obtain 



[PtPt^luud = 2S+ uud (T«,T b ) , [P 3 P?A(% d = \{T*,T b } ds o, 

[^2 + ^i + A (2) ] sW = 2St uud {T\T b ) , [P 3 P?A(% d = \{T\T b } ds {-^N v 

[P^Pi + ^%uu d = 2Sf uud (T",T b ) , [P 3 P?A(% d = ±{T«,T b } ds (aN v + 3), 

[P2 + Pi + ^% uud = 0, [P 3 P?^% d = \{T\T b } ds {N v + 2a). 

(A.21) 

For N v = 3, the function A$ = 2Sj uud {T a ,T b ) is shown explicitly in Eq. (4.2). 



Appendix B. Ultraviolet divergences and 0(p 4 ) operators 

Once we go beyond the order 0(p 2 ) in xPT 1 the number of operators that enter Eq. (2.10) 
increases dramatically. At the order 0(p 4 ), we rewrite the weak Hamiltonian as 

H w = 2V2G F V ud V: s {\ \g 2 70 2 7 + J2 D l 6$]+2\g 8 8 + g' 8 0' 8 + Y,E i 6 ( i ) ])+R.c. , (B.l) 



iVv 



3, (over)complete sets for 



27 ' 



where 0§ , of are the new operators. For TVf 

have been listed in Ref. [60]. The use of partial integration identities makes it possible 
to reduce the number of operators drastically, leading to the lists commonly used in phe- 
nomenology [61]; in our case, however, the use of partial integration identities is not possible, 
since we consider local operator insertions (i.e. 7i w is not integrated over spacetime). 
Generalizing Eq. (2.20), we define the correlation functions now with the LECs added, 



[C 2 7] a \x ,y ) = Jd 3 xJ a 3 y(j a (x)[g 27 O 27 (0) + ^D t 6$(0)]j b (y)) , 
[Cs] a \xo,yo) = Jd 3 xJd 3 y(j a (x)[g 8 Os(0) + g 8 O 8 (0) + Y^E i 6$\0 



The results can be written in the forms 



[C 27 ] a (x ,y ) = A2 b \ g 27 C(x )C(y ) +V 27 (x 0l y ) + S 27 (x , yo) >, 



(B.2) 
(B.3) 

(B.4) 



[C 8 } a \x ,y ) = Afig 8 \C{x )C(yo)+V 8 (x ,y )\ + g' 8 V' 8 {x , yo) + S 8 {x ,y ) j , (B.5) 

where A$ = 2Sj uud (T a ,T b ) in the notation of Eq. (A.21), and Af = {T a ,T b } ds /2. 

The list of operators from Ref. [60] (modulo certain minus-signs) that can contribute to 
C 2 7 at NLO is constituted by the properly projected (cf. Appendix A) versions of: 



a 
a 



(2) 
27 
(4) 
27 



(,'P)ur(.'P)vs j 

(40™. (4.. 



vs ' 



(B.6) 
(B.7) 
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= (^W-(^k-Tr (S) , (B.8) 

{ 2? = i(W^)- r -(Vh- s , (B.9) 

O?7 0) = -(4.)«r(^W w , (B.10) 

^7 1} = -(40 w (^w w ) w , ( B - n ) 

0g 4) = -(W^) W (W^) M> (B.12) 

0? 7 5) = -(W^-CW^- • (B.13) 



Here we utilize the notation 



5 = U X ] + XU* , 7> = i (W - X^ f ) , (B.14) 
C^Udplft , W^ = 2{d^C v + d v £ ll ) , (B.15) 

where % = 2mS/i ?2 = M 2 . As stressed in Ref. [60], not all of these operators are independent, 
however: equations of motion can be used to eliminate 19, 21, and 25, for instance. In the 
following, we keep for generality all the operators. 

The contribution from the 0(p 4 )-constants to Eq. (B.4) reads 



£27(^0, 2/0 ) = 4M 4 P'{x )P'(y ) [D 2 + 2D 19 - 4D 24 - 4D 25 ] + 



+ 4M 6 P(x )P(y ) 



Nr 

D 4 + -j-D 7 - D 20 - D 2 i - F 2 (N f L 4 + L 5 ) 



(B.16) 



Taking into account that C(xq) is finite; that T> 27 (xo,yo) contains the divergences specified 
in Eq. (4.8); that the QCD 0(p 4 ) constants contain the divergence (A = — l/32-7r 2 e) 

N t L 4 + L 5 = N f L r 4 + L r 5 + ^-\, (B.17) 

where L 4 , L\ are finite; and that the 0(p 4 ) constants contain the divergences 

D4 = Dl + g 27 F*\{^±^ 

D 7 = D r 7 +g 27 F 2 \(\ 



D 20 = D r 20 +g 27 F 2 \(- 



D 



21 



D r 2A + g 27 F 2 X^ , (B.18) 



the correlation function C 27 in Eq. (B.4) can be seen to be finite. 

As far as the octet correlation functions are concerned, it is the following types among the 
operators listed in Ref. [60] that contribute to the correlation function Cg at the order we are 
considering: 

0f } = ~(SS) ds , (B.19) 



27 



o 



o 



(2) 
8 

(3) 



6 

O, 



(10) = 
8 

(11) = 
8 

(14) = 
8 

(33) = 
8 

(35) = 
8 

(36) = 
8 — 

(39) = 
8 

(40) _ 



-(S) ds Tr(S) , 
-{VV) ds , 

1*5) £fj,£-iJ,}ds i 

{CftC^dsTr (S) , 

-{£^,d u W^} ds , 
-{C^d^W uu } ds , 



(B.20) 
(B.21) 
(B.22) 
(B.23) 
(B.24) 
(B.25) 
(B.26) 
(B.27) 
(B.28) 
(B.29) 



Again, there are relations between these operators: equations of motion can be used to 
eliminate 33, 36 and 40 [60]. For Nf = N v the contributions from the QCD and weak 
0(p 4 )-constants to Eq. (B.5) read 



Nt 



£s(x ,y ) = 8M 4 P'(x )P'(y ) 
+ 8M 6 P(x )P(y 
+ 4M 6 P(x )P(y )g s F 2 [-7V f L 4 - L 5 ] + 



E\ 2~^ 2 + -^3 + 4i^33 — 4^39 — AE40 



+ 



E\o + 2^11 + ~2^^ 14 ~~ "^ 35 ~~ 2-^36 



+ 



+ 4M e 



dM 2 



P'(x )P'(y ) g' 8 F 2 -7V f L 4 - L 5 + 2(N { L 6 + L 8 ) 



(B.30) 



The results for the divergent parts of Ei can be found in Ref. [60] for iVf = 3 and will be 
given below for general Nf. Taking into account that (in the unquenched case) 

A 



N { L 4 + L 5 - 2(N { L 6 + L 8 ) = N { L r 4 + L\ - 2(N { L r 6 + L r 8 ) + 



4iV f 



(B.31) 



where Lg, Lg are finite, and summing together with the divergences shown in Eqs. (5.9) and 
(5.15), it can be verified that Cs is finite. 



Appendix C. The case N { ^ N v 



For iVf 7^ N v , the set of possible operators is in general larger than for N{ = N v : the only 
restrictions are that the operators be singlets in the full group SU(iVf)^, and have the correct 
transformation properties in the subgroup SXJ(N v )l. At 0(p 2 ) this does not change the 
situation for the 27-plet, but it increases the amount of octets to four in total. Besides Og 
defined by Eqs. (2.13), (2.12), (2.8), viz. 

F 4 



(C^C^) ds + (C,) ds TT{P^) 



(C.l) 
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and 0' 8 defined by Eq. (2.14), there are two additional octets, which we choose to define such 
that they vanish in the limit Nf —> N v : 



Os = ^-[^(l-Pv)^U, 
F 4 



8 = 



(^UTr[(l-^v)^]. 



(C.2) 
(C.3) 



It should also be noted that these operators only contribute starting at the NLO, since at 
tree-level they do not couple to two valence-flavoured mesons. Since for C27 nothing changes 
with respect to Appendix B, we concentrate on the octets in the following. 
The three-point octet correlation function is now of the form 

[C 8 ] ab = [d 3 x[ d 3 y (jg(x) \g 8 8 + g' 8 0' 8 + g 8 6 8 + g 8 d 8 + £ E&f] (0)J$(y)) (C.4) 



= Af lg 8 [C(x )C(y ) +V 8 (x ,yo)\ + g' 8 V' 8 + g 8 V 8 + g 8 V 8 + £ 8 j , 
where V 8 can be found in Eq. (5.5) and V 8 in Eq. (5.14). The new functions read 
V 8 (x ,y ) 

£•8 (^0,2/0) 
where we have defined 

l A (x ,yo) = G(0;M 2 )P'(x )P'(y )- 



F 2 M 2 



1 F 2 M 2 
■-T>27(x , y ) H — 1a{xq, yo) 



N v l B (x ,y ) - 1a(x , yo) 



M 2 r 

— ^-P(x - Vo) [B(xo) + B(y ) 
+ ^P'{x -y )[B'(x )-B , (y ) 



+ 
+ 



+M 4 f dr B(t)P(t - x )P(t - y ) , 
J 



l B (x ,yo) = E(0;M 2 )P'(xo)P'(yo)- 

M 2 [~ 
— ^-P(x - yo) [B{xq) + B(y ) 



+ 



+ 2^0 - yo) [B'(x ) + B (x ) - B'(y ) - B (y ) 



+ 



+M 2 [ dr 
Jo 



(C.5) 

(C.6) 
(C.7) 



(C.f 



M 2 B(t) + -Boo(r) P(t - x )P(r - y ) , (C.9) 



and the notation follows that in Eq. (5.5). The divergent parts read (in the unquenched case) 



V 8 (x ,y ) = V r 8 (x ,yo) + ^(N[-N v ) M«P(x )P(y 



(CIO) 



^(^o, yo) 



^sO^o, yo) + 



F 2 \ 



(l - ^) [M e P(x )P(yo) ~ M i P'(x )P'(y ) (C.l I ) 
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where Pg, Pg are finite. 

The list of operators contributing to £g for Nf / N v is also much longer. We will not provide 
any systematic classification of all the possibilities, but only list the additional operators that 
are needed for cancelling the ultraviolet divergences at NLO. Using the same notation as in 
Appendix A [Tr v (...) = Tr (P v ...)], we need 

(C.12) 
(C.13) 
(C.14) 
(C.15) 
(C.16) 

(C.17) 
(C.18) 
(C.19) 
(C.20) 



u 8 — 


-(SP v S) ds , 


o (2 ' ] = 


-(S) ds Tr v (S) , 


^(10') _ 
u 8 


{<5 ', £■ fj,Pv £ ^} ds j 


rfW") - 
u 8 — 




u 8 — 


2 (Ai{-fv> S}£/i)ds i 




{CpPyCuUTr (S) , 


^8 — 


(CpCfJdsTry (S) , 


/n( 35 ') - 
^8 — 


-(c^d.w^ + dvW^p^ 


r,(39') _ 

u 8 = 


-(W^P v W^) ds . 



With these definitions, we get: 



N v 



£s(xo,y ) = 8M 4 P'(x )P'(y ) -E l - E v - -^E 2 - -^E 2 , + E 3 + 



2 z 2 
+ 4E 33 -4(E 39 + E 39 ,)-4E 40 
1 



+ 



+ 8M b P(x )P(y ) [E 10 + E w , + E w „ + -(E u + E 1V ) + 

N N v 

+ ~Y^ E U + E 14') + ~Y E W ~ 2 ( E 35 + ^35') ~ 2E 36 

+ 4M 6 P(x )P(y )g 8 F 2 [-JV f L 4 - L 5 ] + 



+ 



+ 4AP 



dM 2 



P'(x )P'(y ) g' 8 F 2 -NfU - L 5 + 2(N { L 6 + L s ) 



(C.21) 



Like for C27, the part C(xo)C(yo) in Eq. (C.5) is finite, while the other parts contain 
divergences. More precisely, Vg, P 8 , Pg, Pg are of the forms shown in Eqs. (5.9), (5.15), 
(C.10), (C.ll), the combination NfL^ + L5 of the form in Eq. (B.17), while the combination 
on the last line of Eq. (C.21) is of the form in Eq. (B.31). Moreover, writing 



Ei = El + 



F 2 X 



(gsm + aWi + hm + Mi) 



(C.22) 



where E\ are finite, the coefficients r/j, r/-, fji, fji can be derived with the method of Ref. [60]; 
they are listed in Table 1. Summing together, all the divergences cancel in Cg, as they should. 
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1' 



10 



Vi Vi Vi 



i ^v v -rz^ 16 AN{ j 4 16 4A r f 8 -r 4JVf 



1 L. n 1 

8 4JV f u 47V f 



2 f7V +21^7 -- ^7 - ^7 

^ 8 4JV f u 8 4AT f 

3 



yiy v -t-z,jy 32 -r leNf J 8 32 16 ^± ^ y 



10' J- + ^£ o 

iu 16 ~ 16 u 16 16 

in" _A — 

16 ' u 16 

Nt 37Vv 

8 8JV f 



11 (JV v + 2)(-g^) 

11' 3 , JVjf n _JVf _3 

±x 8 8 u 8 8 

14 i 

14' I -\ 

14" -A ^ 

33 (^ + 2)^-3^) ^ -^(l-f 



35 (iV v + 2)(-^) ^ B 
35' ^ 

36 

39 (JV v + 2)(-^) & 
39' ^ —t^ 

40 



Table 1: The coefficients that appear in Eq. (C.22), in the unquenched case. 

The results for E^ that need to be used for the case Nf = N v can be obtained from Table 1 
by summing together the coefficients with the same "numerical" index that then correspond 
to the coefficients of the operators in Eqs. (B.19)-(B.29): (E\ + £?i')jVf=JV v f° r Eq. (B-19), 
(E2 + E2')n { =n v for Eq. (B.20), etc. It can immediately be seen that the divergent parts 
proportional to gs and gs cancel in these sums, as has to be the case. 

We finally comment on the quenched limit, corresponding formally to Nf — > but iV v 
fixed. We have seen that for Nf 7^ -/V v additional operators in general appear, as elaborated 
in Ref. [32]. However, it is easy to see that the functions 2a, 2b that appear in Eqs. (C.6), 
(C.7), vanish in the e-regime. Therefore the coefficient gs does not contribute in Eq. (C.5) in 
the e-regime. Moreover, T>s is determined by the same function D27 that appears in T>$ (cf. 
Eqs. (5.5), (C.6)). In particular, the normalised three-point function defined in analogy with 
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Eq. (5.10) obtains for Nf — > the form 



A 



ab 



98 ~ (98 - 9a) Tpj\2 (p »/?i-P*0o) 



+ f8[ C 8]norm( 2; 0,yo) • 



(C.23) 



We observe that quenched functional behaviour only appears in the part [Cg]^ rm (a;o, yo) (cf. 
Eq. (5.23)), and can thus be eliminated by disentangling the contributions related to g' 8 , 
just like in Sec. 5.2. Moreover, it can be verified that at the NLO in the e-regime, the 
coefficients gs, gs do not contribute to the correlation function considered in Sec. 5.3, such 
that g' 8 can be separately determined just like there. The remaining terms in Eq. (C.23) can 
be disentangled in principle by monitoring the volume dependence, from which it should be 
possible to determine the "physical" coefficient g%. 



Appendix D. Correlation functions for N v = 4 

It has been argued recently that many of the mysteries related to the AI = 1/2 rule can be 
studied particularly cleanly [both from the conceptual and from the practical point of view] 
by considering the SU(4) symmetric situation, i.e. N v = Nf = 4 [20]. We discuss here how 
our predictions can be converted to apply to that situation. 

Rather than 27, 8, the dimensions of the relevant irreducible representations are 84, 20 for 
N v = 4. The corresponding operators are obtained like the 27 for N v = 3, but by using the 
projection operators P%P\ in Eq. (A. 4), with a = +1 for the 84 and a = — 1 for the 20. 
Following the notation in Ref. [20], the corresponding operators are denoted by [Oi^mr 

The three-point correlation function we are interested in now takes the form 

[CiC>o,2/o) = jd 3 x |d 3 y(^(*){^[0^^ (D.l) 

The 0(p 4 ) weak operators Oi here have the same chiral structures as the 27-plets of Ref. [60], 
listed in Eqs. (B.6)-(B.13) of Appendix B, but each of them comes in two variants after the 
valence flavour projection, corresponding to a = ±. The result can be written in the form 

[CiCuvixo, yo) = jtf [c(x )C(y ) + aV 27 (x , y )] + S°{x , y )} , (D.2) 

where = ^ suv (T a ,T b ), with the function S? suv (T a ,T b ) given in Eq. (A. 14). The 

function T>27{xo,yo) is identical to the one for the 27-plet in Eq. (4.6). 

The functions <S°"(xo,yo) in Eq. (D.2) contain the contributions of the 0(p 4 ) low-energy 
constants, beyond those already contained in the factorized term C(xo)C(yo): 



£ a (x ,y ) = 4M 4 P'(x )P'(y ) [D% + 2D^ 9 - 4D% 4 - AD^} + 



+ AM 6 P(x )P(y ) 



Dl + - D5o - Dfi - F 2 (N { L A + L 5 ) 



(D.3) 
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Taking into account that C(xq) is finite, that the V27 contains the divergences in Eq. (4.8), 
that the QCD 0(p 4 ) constants contain the divergence in Eq. (B.17), and that the weak 0(p 4 ) 
constants contain the divergences 



Dl = Dl r + <7?F 2 A( 




(D.4) 



D7 = D^ + gfF 2 \Q, 
Dl, = D%+g<[F 2 \^), 



(D.7) 



(D.5) 



(D.6) 



the correlation function C± in Eq. (D.2) can be seen to be finite. 
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